Universal quantum control in irreducible state-space sectors: 
application to bosonic and spin-boson systems 
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We analyze the dynamical-algebraic approach to universal quantum control introduced in P. Zanardi, S. Lloyd, 
quant-ph/0305013 The quantum state-space Ti encoding information decomposes into irreducible sectors and 
subsystems associated to the group of available evolutions. If this group coincides with the unitary part of the 
group-algebra CIC of some group K, then universal control is achievable over the /C-irreducible components 
of H. This general strategy is applied to different kind of bosonic systems. We first consider massive bosons 
in a double-well and show how to achieve universal control over all finite-dimensional Fock sectors. We then 
discuss a multi-mode massless case giving the conditions for generating the whole infinite-dimensional multi- 
mode Heisenberg-Weyl enveloping-algebra. Finally we show how to use an auxiliary bosonic mode coupled to 
finite-dimensional systems to generate high-order non-linearities needed for universal control. 
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I. INTRODUCTION 

Over the last few years quantum information science taught 
us how to take advantage of information encoded in the states 
of a quantum system 1 1]. A fundamental requirement is then 
the ability to process such information in the most general 
fashion by resorting to the physically available interactions. 
When arbitray transformations (or approximations of) over 
the system state space H are dynamically realizable one says 
that universality is achieved. To this end one typically in- 
troduces a set of elementary building blocks, e.g., quantum 
logic gates, whose combinations generate the whole set of 
desired dynamical evolutions. A very important example is 
provided by quantum computation |1]. There the prototype 
information-encoding system is provided by a collection of 
two-level systems i.e., qubits, such that each of those subsys- 
tems is completely controllable i.e., any SU(2) rotation is re- 
alizable, and an entangling two-qubit transformation is avail- 
able H,H,li. 

Unfortunately there are several physically relevant situa- 
tions e.g., spin-based quantum computing, in which such a 
goal is extremely difficult to achieve practically. In those cases 
the lack of easily available dynamical resources forces one to 
consider more sophisticated ways of encoding and manipulat- 
ing information. For example it might be the case that inde- 
pendent single qubit-operations are very difficult to realize in 
view of the small spatial separation between qubits. In this 
case just global operations affecting the quantum register as 
a whole are readily available; these latter clearly do not gen- 
erate the full set of transformation over H. Another class of 
examples of this state of affairs is provided when all naturally 
available interactions commute with some observable e.g., to- 
tal spin. It follows that the state-space splits in mutually in- 
accesible sectors corresponding to the different eigenvalues of 
such an observable. 

This kind of problem led several authors to develop the con- 
cept of encoded universality O Is III l2l Gal • That was mostly 
done by explicitly working out a few important and non-trivial 
istances of a such a notion. Roughly speaking the idea is that 



even in presence of a limited set of available control resources 
one can find out a in invariant subspace C of H over which 
universality is achievable. More recently a general dynamical- 
algebraic framework has been suggested to underlie all these 
former examples and to provide room for further generaliza- 
tions and applications 1 1 (j . 

The key feature of the techniques advocated in Ref. fioll is a 
massive use of tools borrowed from group representation the- 
ory Hill . Assuming that the set evolutions one can enact forms 
a group i.e., the control group Ua, the state-space splits ac- 
cording the associated irreducible representations (irreps). In 
such invariant subspaces Ua typically still does not allow for 
universal control, nevertheless that goal can be accomplished 
when one has at disposal generic Hamiltonians belonging to 
the group algebra of Ua- It follows that if one recognizes 
that the controllable interactions belong to the group-agebra 
of some hidden non-abelian group /C universal controllability 
is realizable over the whole family of subspaces and subsys- 
tems associated with the irreps of K, I Kill . In this approach an 
important role is played by the use of symmetries that allow 
one to build dualities between different kind of encodings and 
manipulation strategies. This crucial role is also at the basis 
of the deep conceptual connection between universal control- 
lability strategy proposed in 1 1 ] and the theory of quantum 
noise correction and avoidance ll^ . fPill . ffl . 

In this paper we shall further develop and apply the general 
lines of Ref. llOll . The relevant conceptual and mathematical 
framework will be illustrated in a fully detailed manner and 
novel physical examples worked out. A major emphasis will 
be given to bosonic systems, both massive ones (with associ- 
ated finite dimensional irreducible sectors) and massless one 
(infinite dimensional Fock space). 

The paper is organized as follows. In section II will be in- 
troduced the basic terminology and quantum-control notions. 
In section III we will lay down the general group-theoretic 
formalism at the basis of our approach. In section IV and V 
we will analyze the controllability of a pair of coupled mas- 
sive bosonic modes in terms of the associated SU (2) struc- 
ture. In section VI the massless bosonic case with infinite- 
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dimensional state space will be considered. In section. VII we 
will discuss an auxiliary bosonic mode technique for generat- 
ing highly-non linear interactions. Finally section VIII will 
contain the conclusions. 



II. PRELIMINARIES 

The starting point is the set of physically available interac- 
tions that act over the quantum state-space H: 

1 A ■= {H{\)} XeM C End(H) (1) 

where M. is the set of possible values of the control param- 
eters A. The physically realizable control processes can be 
described by A(t) £ Va, i-e. by a set of .M-valued functions 
(paths) that represent the evolutions of the control parameters 
that are experimentally achievable. Thus the unitary quantum 
evolutions one can enact are those that are governed by the 
time-dependent Hamiltonians H(X(t)): 

U(X) = Texp(-i [ H(\(t))dt) (\eV A )- (2) 

J JR. 

The pair (Ta,Va) describes the physical resources available 
in the given experimental situation. 

In the following we will assume that if U is an allowed 
evolution, then W is allowed as well; we also assume that 
the trivial evolution, i.e. U = 1, is an allowed one. With 
this assumption the set of unitary transformations that one can 
generate in the experimental situation described by (Ta, V a) 
has the structure of a group Ua ■ 

We now give a brief description of the properties of the 
group Ua that depend on the paths Va that can be real- 
ized. We start by describing a well-know result of quantum 
control theroy. Suppose that for any two arbitrary operators 
A, B 6 Ta one is allowed to apply sequences of pulses 
cxp(±iA5t), exp(±iBSt), than the following relations hold 

e±i ASt e ±iB5t = et (±A±B)8t + 0(Jt 2j. (3) 
e iASt e iBSt e -iASt e -iBSt _ e i(AB-BA)St 2 + Q(St 3 ). (4) 

This means that if one can drive the control parame- 
ters along arbitrary paths in M., then one can generate, 
by iteratively applying the sequences and @, all the 
possible linear combinations of Hamiltonians of the form 
±iA,±i[A,B],±i[[A,B],C], ect. , where A, B, C, ... € I A 
0]. The result can be stated in terms of unitary operators as: 

U A = e CA , (5) 

where La is the set of Hamiltonian generated from the set Ta 
by commutation. La it has the structure of Lie algebra. On 
the contrary, if we suppose that there are restrictions on the 
experimentally realizable paths Va then the allowed unitary 
transformations are such that Ua C c Ca . 

From the point of view of the control theory the question 
that naturally arises is the following: given the pair (Ta, Va), 



what kind of control one has on the physical system S under 
study? One has universal control over S when it is possible to 
enact all the unitary operators U(TL) where H is the quantum 
state-space of S. Now, depending on the pair (Ta,Va) there 
can be different situations: 

• if Ua = U (fi) universal control over H is achieved; 

• if Ua is dense in U(H) then universal control over H 
is still achieved since it is possible to simulate any uni- 
tary evolution with arbitrary accuracy by means of the 
available resources 

• if Uai is not dense in U(fi), then it is not possible to 
have universal control over H; in this case one can in- 
dtroduce the notion of encoded universality. 

When the pair (Ta,Va) does not allow for universal control 
over the whole state-space H it is natural to find out whether 
there are subspaces Ci C H such that: 

U A \ Ci =U{Ci). (6) 

In this case we say that we have encoded universality, or 
that Ua is (^-universal. The subspaces C;'s are referred to as 
codes, in the framework of quantum information theory these 
are the subspaces in which one can encode the information to 
be processed. 

It is interesting to note that the C,*'s need not to be invari- 
ant subspaces for Ta, in fact the encoded universality can be 
achieved even by allowing the state of the system to temporar- 
ily leave the code subspace during the time evolution. A sim- 
ple example of this situation is given by the subspaces of a 
code C: if Ua is C-universal than it is universal over all the 
subspaces C' C C. In this case the evolution starts and ends 
in a given subspace, althought other intermediate states (an- 
cillce) that belong to (C) might be reached during the time 
evolution of the system. 

It is clear that, from the point of view of both quantum con- 
trol theory and quantum computation theory, it is fundamental 
to find a general way to define which kind of control one has 
over a given quantum system. In the following section we give 
a general procedure to address this problem. 

III. GENERAL FRAMEWORK 

The idea is to find out under which conditions the physi- 
cally available interactions allow to perform universal control 
over the system S or, if this is not possible, to give a proce- 
dure to identify the codes over which one can have encoded 
universality. We start by giving a concise description of the 
mathematical ingredients involved and at the same time giv- 
ing a first physical example of their realization. 

The example is given by a three-levels system which is well 
known in atom-lasers physics . The Hilbert space of the sys- 
tem is given by span{|l), |2), |3), |e)} where the first three 
states are the degenerate (atomic) ground states, |e) is the first 
exited state; AE is the energy difference between the excited 
and the ground state. The system is coupled with controllable 
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external potentials which are lasers characterized by the Rabi 
frequencies fii, 51 2 , ^3 G C. Since the lasers are detuned with 
respect to the energy gap AE, the ground states are coupled 
to one another by a second order process, but they are never 
coupled to the excited state. The Hamiltonian of the system 
can be written as 

H = Qi\l){2\ + n 2 |2)<3| + n a |3)(l| + h.c. (7) 

and the evolution of the system takes place in the degenerate 
ground state subspace, i.e. H = span{|l), |2), |3)}. 

We now start by giving the basic ingredients of our control 
picture: 

• the system S, whith state-space H, and the correspond- 
ing group of unitary operators U(TL); 

• the set of physically available and controllable interac- 
tions 1 A := {£\ \Hi}; 

• the Lie algebra Ca generated from 2^4 := XiPn}; 
by commutation; 

• a discrete (continuos) group JC and its (irreducible) rep- 
resentations; 

• the group algebra C/C generated by JC; 

• the unitary part of the group algebra UCK,. 

In our first example we have already seen that Tt = 
span{|l), |2), |3)}, thus the unitaries over H belong to U(3). 
Since we suppose to be able to control the three lasers sepa- 
rately, the set of avalable interactions is given by: 

I A := {fii|l)(2| +h.c.,Q 2 \2)(3\ + h.c.,n 3 \3)(l\ + /i.e.}. 

Suppose now that JC is a discrete group. We start by giving a 
description of the group algebra C/C |22],[23]. This is defined 

as 

C/C = {x I x = J2 X a9 ; A s G C}, (8) 

that is, it is the vector space of all linear complex combina- 
tions of elements of JC. Its dimension is dim^ = |/C|, i.e. 
the number of elements of JC. In order to complete the defini- 
tion of the group algebra one has to define: i) a binary oper- 
ation (product) on it, ii) a unary operation (conjugation). The 
product law can naturally be defined by means of the group 
multiplication law; if x, y G C/C, that is x = X^eje X a9 anc ^ 
V = Ehe/c X hh then: 

xy= X g X h gh 

g,h£K 

belongs to C/C since gh G /C, Vg, h . The conjugation is sim- 
ply defined Vx G C/C as: 

= E ^ 

g&K. 



where A g is the complex conjugate of X g , and g^ 1 is the in- 
verse of g with respect to the group law. 

We will be interested in the representation of the elements 
of JC as operators over a Hilbert space Ti, i.e. p : JC — > C(TL) 
a linear map onto the group of linear operators on TL. It is 
interesting to notice that by means of the representation p one 
can give a representation of C/C over TL. In fact Vx G C/C 
the operator p(x) — J2 g eic X aP(d) belongs to C{TL). Thus, 
if H has finite dimension, from now on we can think of the 
elements of JC and C/C as complex n x n complex matrices, 
where n = dim7i. 

In our example we can choose the group JC = V3, i.e. the 
group of rigid symmetries of a triangle. We have that 

D 3 = {1,R, P, RP, PR, RPR}, (9) 

that is D3 has six elements and it can be written in terms of the 
two generators R and P , which represent a rotation of 27r/3 
and a reflection, respectively. If we now represent the vertices 
of the triangle with the basis states |1), |2), |3) we see that we 
can give a three dimensional representation p of D 3 by writing 
R and P in terms of operators over H = span{ 1 1) , 1 2) , 1 3) } in 
the following way: 

p(i?) = |l)(2| + |2)(3| + |3)<l|; p(P) = |l)(2| + |2)(l|. (10) 

By means of the same representation we can give a represen- 
tation of C2?3, so that we can think about the elements of D3 
and C2?3 as 3 x 3 complex matrices. 

It is now interesting to have a closer look to some of the 
elements that compose C/C. In order to do this we first notice 
a couple of nearly obvious facts: i) if x, y G C/C then also 
xy—yx belongs to C/C, i.e. the commutator of two elements of 
C/C is again an element of the group algebra; ii) if we consider 
unitary representations of JC, p : JC — > U(H), as it is in 
the case of[H)| then the representative of the conjugate of an 
element of C/C corresponds to the Hermitean conjugate of the 
representative of the element itself i.e., p(x* ) = p(xy. 

The first subset of C/C that is fundamental for universal con- 
trol is the unitary part of C/C defined as: 

tt(CAC) = {x G C/C \xx* = x*x = 1}; 

that is it is the set of elements of C/C that are represented by 
p as unitary operators over Ti. It is easy to verify, using the 
multiplication law of JC, that this set is actually a group. Since 
we suppose p to be a unitary representation of JC we have 
p(IC) C p(U(<CJC)). The other important subset of C/C is the 
set of all the anti-hermitian operators: 

u(CJC) = {x e CJC \ x = ~x*} 

this is actually a Lie algebra, furthermore we have that 
U(CJC) = exp(u(C/C)), that is the unitary elements of the 
group algebra can be built out of the elements of it (C/C) via 
exponentiation. In our example we have that 

u(€V 3 ) = iaP - f3R + f3R*, (11) 

with a G H, j3 G C 
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We now give some standard results of the theory of groups 
representation that will be useful in the following paragraphs. 
Suppose p is a representation of the group JC on Ti. Then p 
is a reducible representation if there is at least one subspace 
Ti' C Ti that is invariant with respect to the action of the 
elments of JC, that is Ti 1 is such that V x E JC and V \4>) € Ti' 
we have that p(x)\4>) e Ti.'. In this case the restriction p\ w 
is itself a representation of JC on Ti' . If the only invariant 
subspaces of Ti with respect to p are Ti itself and {0} then the 
representation is called irreducible. 

In general, if p is reducible it induces on Ti the following 
decomposition: 

U = (BjeiHj, (12) 

that is Ti can be written as the direct sum of the orhotogonal 
invariant subspaces Tij, on which the group acts irreducibly 
(X is the set of indexes that label the various subspaces). Thus 
the restriction pyn , of the representation to any of these sub- 
spaces is an irreducible representation of JC that has dimension 
dj = dimTij. Notice that in Eq. d!2t the sum is over possi- 
bly equivalent irreps. If the label J is meant to represent just 
inequivalent irreps, then one has to introduce factors <D™ 7 tak- 
ing into account the multiplicity nj with which the J-th irrep 
appears 1 10[. 

In the case of the three-level system, it can be easily 
checked that, given the representation jlOl l of D 3 , the Hilbert 
spaced = span{|l), |2), |3)}canbe written asTi — Tii®Ti 2 , 
where 

Wi=span{(|l) + |2) + |3))/V3}, (13) 

3 3 

Ti 2 = span{^]exp(i27rj/3)|j),^exp(i47rj/3)|j)}. 

j=i j=i 

(14) 

Thus, Tii and H2 are invariant subspaces, and we have that the 
restriction p\ ni (p\n 2 ) °f 03 to Hi (H2) gives a one (two) 
dimensional representation of D 3 . 

We can now give the result of representation theory that is 
fundamental for the universal control. Suppose p is an irre- 
ducible representation of JC over Ti, then, with respect to the 
same representation, 

U{€K)=U{H). (15) 

This means that if, with the available and controllable inter- 
actions I a, we are able to generate by commutation (see (0 
and @) all the hermitian elements of u(CJC) (Hamiltonians), 
then, since U(€JC) = exp(u(<CJC))) — U(Ti), we can have 
the universal control over the system S. 

The result can be extended to the case in which p is a re- 
ducible representation of JC over Ti. In this case for any ir- 
reducible subspace Tij of the decomposition d!2i . we have 
that 

U(CJC) lHj =U{Hj). (16) 

This means that if, with the available and controllable inter- 
actions Ia , we are able to generate by commutation (see 



and @) all the anti-hermitian elements of w(C/C), then we 
can generate all the unitary operators over each Tij and so 
we can have encoded universality over each of the irreducible 
subspaces. 

In the case of the three-level system we have that H2 is 
an invariant subspace for the representation J10I of D 3 , and 
the restriction p\n 2 of (II Oi to H2 is irreducible. This means 
that, since with the available interactions we are able to gener- 
ate the whole u(€V 3 ), see {TT}, we have that U(€V 3 )\ H2 = 
U(Ti2), i.e. we have encoded universality over the code 
C2 = H2 (and the same is obviously true for the one dimen- 
sional invariant subspace Tii). 

We are now in the position to summarize our conceptual 
framework for universal quantum control. 

The basic ingredient is the set of physical interactions Ta 
available and controllable for a given system S. In order to 
check if the controllable interactions allow to have universal 
control over Ti, one has to: 

• find a group JC and a representation p over Ti; 

• find the decomposition Ti = (BjTij with respect to p; 

• check if the set of available interactions Ta allow, via 
commutation (see and Q) to generate the whole Lie 
algebra u{€JC). 

Then we have two cases: i) p is irreducible: then one has uni- 
versal control over the whole state-space Ti; ii) p is reducible: 
then one can achieve encoded universality over all the irre- 
ducible subspaces Tij (codes) of the decomposition (I12> . 
A first important example for the group JC is 

JC=U A ^eMC A ). (17) 

In order to find out if it is possible to have encoded universal- 
ity, one has to look at the irreducible representations of Ua- 
In this case one thing should be clear: although Ua acts irre- 
ducibily on each of the subspaces Tij of the decomposition 
O l this does not imply that U(VUa)\h , = Utfij). In fact, 
to have encoded universality in the irreducible subspaces one 
has to prove that 

L A =u{<CU A ), (18) 

i.e. the Lie algebra generated from the set Ia via commuta- 
tion coincides with the Lie algebra of all the anti-hermitian 
elements belonging to the group algebra U (CUa)- 

Before concluding this section it is worthwhile to stress that 
the dynamical-algebraic machinery we developed allows one 
to sistematically build classes of (to begin with) formal exam- 
ples and then to look whether they admit some natural physi- 
cally realization. The way to proceed in this case is: 

i) choose a group JC and a representation p over an Hilbert 
space Ti; 

ii) find the decomposition Ti = QjTij with respect to p; 

iii) find the minimal set of interactions I a that allow via com- 
muation (see and Q), to generate the whole Lie algebra 
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iv) find a physical system whose state space is isomorphic to 
H and for which one has the control of the minimal set of 
interactions I a- 

Then, for the given physical system, one can generate all the 
unitary evolutions for each of the subspaces Hj. 

In the next sections we will analyze some physical mod- 
els and show how the general theory applies to them. In 
Ref. IK ill the emphasis was mostly on spin 1/2 systems i.e., 
qubits, with reducible action of the available evolutions group 
Ua- Here we will be mainly concerned with different kind of 
bosonic systems relevant to quantum information processing. 
For those systems the group of available evolutions will have 
an irreducible action over the physical state-space. Accord- 
ingly those examples cannot be regarded, strictly speaking, as 
instances of encoded universality. 



IV. BOSONIC DOUBLE- WELL 

In this section we focus on a well known and well studied 
model: the Bose-Hubbard (BH) model for two boson wells 
filled with a fixed number of bosonic particles N 1 17]. In this 
model, the two wells are given by an external potential ap- 
propriately shaped. The particles occupy the lowest energy 
levels in each site and can be allowed, by lowering the bar- 
rier between the wells, to tunnel from one site to the other. 
We suppose that the offset of the ground state of a site with 
respect to the other can be controlled by applying an addi- 
tional external potential. Furthermore, the particles can inter- 
act among each other when they are in the same site via a two- 
body scattering process (self-interaction). The BH model has 
been recently used to investigate many interesting properties 
of different systems going from arrays of BEC where the sites 
form an infinite (periodic) lattice 111 811 . to systems in which the 
number of sites is finite (dimer, trimer,graph structures) fiH ). 
Moreover, this model has recently been used to describe new 
possible schemes for quantum computing 12011 . According to 
this model the Hamiltonian of the system composed of two 
sites filled with N particles can be written as: 

(19) 

where: the indexes 1,2 label the sites; c\ and a ([ci,c|] — 
1) are bosonic creation and annihilation operators that cre- 
ate/annihilate particles in the local well i; rii = c\ci is the 
corresponding occupation number operator. The parameters 
7i can be used to model the offsets in the ground state en- 
ergies of the different sites.. The parameter e controls the 
nonlinear two-body interaction between the particles in each 
site, while the parameter r control the tunneling processes be- 
tween the sites. In general the system lives in the Fock space 
Hf = h® L , where h = {|n)}„ = o....oc is the state space of a 
single quantum harmonic oscillator and L — 2 is the number 
of bosonic wells (sites). If the number of particles N is fixed, 
as we shall assume, the corresponding Hilbert space is 

H N = span{|N - n) 1 \n) 2 } n =o,..,N, (20) 



where \N — n)i|n)2 is the basis state corresponding to the 
situation in which the are N — n bosons in the first well and 
n bosons in the second one. Hn is a subspace of Hf with 

dimension djv,L := y ^L^l ^J' 

In order to study the control properties of the system it is 
useful to resort to the Schwinger picture of spin operators 1 2 1 ] , 
by means of which it is possible to give the following two- 
bosons realization of the spin operators X, Y, Z: 

X = c\c 2 + 4 c i; Y = £(4ci - c\c 2 ); Z = 1 2 . 

These operators are the generators of the Lie algebra su(2) 
and satisfy the commutation relation [^4, B] = itABcC, 
where A,B,C£ {X, Y, Z} and cabc is the totally antysim- 
metric tensor. In this representation the Hamiltonian dl9l can 
be rewritten as: 

H = lA Z + tX + 2sZ 2 + ls N + eN 2 /2, (21) 

where7 j4 = (71-72), 75 = (71 +72 - 2e)/2. N = n 1 +n 2 
is the total number of particles operator and, as long as the 
number of particles in the wells is fixed, it is proportional to 
the identity operator. 

If we now suppose to be able to control the parameters 7^, 
the offset between the ground states of the two wells, r, the 
tunneling rate of the particles and e, the value of the self- 
interaction among the particles in the wells, we see from (12 1 1 
that our set of available interactions is Ia = {X, Z, Z 2 }. 

With this picture in mind, the group K. that naturally arises 
is SU{2) — exp[sw(2)]. The irreducible representations of 
this group are labeled by the index J £ {0, 1/2, 1, 3/2, 2..} 
and the corresponding operators act on Hilbert spaces Hj of 
dimension dj = 2 J + 1. If we now suppose that the num- 
ber of particles in our sistem is fixed to N we see that the 
corresponding Hilbert space Hn is isomorphic to Hj with 
J = (d N ,2 - l)/2. 

Since SU(2) is a Lie group the Lie algebra corresponding 
to the unitary part of the group algebra CSU(2) is the uni- 
versal enveloping algebra of su(2). The latter is the vector 
space spanned by {A a B b C c } a ,b,c=i,..,oo {A ^ B ^ C £ 
{X, Y, Z}), i.e., is the vector space of the polynomials of any 
order in the su(2) generators X, Y, Z endowed with usual Lie 
bracket. Thus, in order to prove that the set of available in- 
teractions allow us to have universal control over each of the 
Hat's, we have to prove that by controlling the set of interac- 
tions Ia = {X, Z, Z 2 } we can generate, by commutation the 
whole universal enveloping algebra of su(2). It turns out that 
this is true and it will be proved in the next section. We can 
then say that we can generate the whole set of unitary opera- 
tors U{H N ) O- 

Once that the possibility to achieve universal control over 
the Hat's has been proven, one has to give for the desired 
value of N an explicit representation of the operators that 
form a basis for the J = (g?at,2 — l)/2 = N/2 representa- 
tion of su(2). 

We work out explicitly the simplest example that corre- 
spond to N = 2. In this case we have that Hn=2 has dimen- 
sion c?2,2 = dj—i — 3. The whole group of unitary operators 
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is then U(3) = exp(it(3)) and the canonical basis for u(3) is 
given by the following 9 elements: 

X nm = \n){m\ + \m){n\, n < m 
Y nm = i(\n){m\ - \m){n\), n<m 
E nn = \n)(n\. 

(22) 

where n, m 6 {0, 1, 2}. In practice, in order to prove that is 
possible to generate all the unitary operators it sufficient to 
find a basis for the Lie algebra sit (3) which is the algebra 
of the 3 • 3 complex anti-hermitian matrices with trace zero, 
whose dimension is 3 2 — 1 Hal . full . Below we write a pos- 
sible set of elements of the enveloping algebra of sit (2) that 
form this basis, each expressed in terms of the basis elements 
of u(3) 





X 


= \/2(Xoi + Xl 2 ) 




Y 


= V2(Y 01 + Y 12 ) 




Z 


= 2(E 2 2 — E no ) 


X 2 


-Y 2 


= 4X 02 


X 2 


-Z 2 


= 2(Xoi + 2Eu — Eoo — E22 


XY ^ 


YX 


= 4F 02 


ZX ^ 


-XZ 


= 2V2(X l2 -X m ) 


ZY - 


VYZ 


= 2V2{Y 12 -Y 01 ), 



(23) 

These relations show explicitly that by controlling elements 
of the eneveloping algebra of the bosonic su(2) up to degree 
two, universality is achieved in the two-boson sector. 

V. ENVELOPING ALGEBRA OF SU(2) 

This section has a prevailing technical nature. The reader 
more interested in following the conceptual stream of the pa- 
per may want to skip it at a first reading. The goal here is to 
prove that it is possible to generate all the enveloping algebra 
of su(2) by starting form the given controllable interactions, 
i.e. Ta = {X, Z, Z 2 }. The procedure is the following: i) 
generate all the elements of order < 2; ii) generate the ele- 
ments of higher order by using a standard procedure that al- 
lows, given all the elements of the enveloping algebra of stt(2) 
of order x > to generate all the elements of order x + !■ We 
first describe how to obtain the elements of order x < 2. Ob- 
viously Y = —i[Z, X], The term ZY is obtained from 

[Z 2 , X] = i(ZY + YZ) = 2iZY - X. (24) 



This first example allows us to describe an obvious an general 
rule that will be used extensively in the general procedure de- 
scribed below. In general the result of the commutation is not 
in the desired form, for example the elements X, Y, Z com- 
pare with the required power but they are not in the required 
order. The ordering process requires the use of the commu- 
tation relations of su(2); these allow to write the starting ele- 
ment, ZY in[24l in terms of a new element of the same order 
in which the position of two factors has changed, and of a 
new term of lower order (X). By subtracting the lower order 
term from the result of the commutator we get the desired one 
(YZ). In the same way we can write the terms ZX and XZ. 
The next step is to obtain the terms X 2 and Y 2 . This can be 
done by respectively adding and subtracting iZ 2 to the result 
of the commutators: 

[Y,ZX} = -i(Z 2 - X 2 ); [X, ZY] = i(Z 2 — Y 2 ). (25) 

The remaining terms XY and YX can be obtained, for ex- 
ample, by applying to [X 2 , Z] the same procedure of ordering 
and subtraction seen for J24i . 

We now focus on the element of order > 2. In order to de- 
scribe the general procedure we use the operators A,B,C <E 
{X, Y, Z} such that A ^ B ^ C. A first obvious basic rule 
that can be applied to generate new elements of the envelop- 
ing algebra is to properly commute elements already available 
at a certain step of the procedure. In general, the result of 
the commutation will be a linear combination of the desired 
new element of order x + 1 with elements of the same or of 
a lower order. Thus, the new element will be obtained simply 
by applying an ordering process on terms of order x + 1 an d 
by subtracting the resulting lower order terms. Another im- 
portant part of the procedure will entail the use of the Casimir 
relation that is satisfied by the elements of su(2): 

J{J+ 1) = A 2 + B 2 + C 2 , (26) 

where J is the label of the representation. This relation will 
be used to build terms of order x + 1 out of terms of the same 
order that are already available. 

A usful relation that will be extensively used is the follow- 
ing; let O be a generic operator, then: 

[A k ,0}= A p [A,0}A q . (27) 

p,9>0,p+q— k— 1 

We now describe the first step of the procedure with which 
we can obtain elements of the type A x+1 . We first write: 



[AX^B.CA] - [A X ~ X C, BA] = is BC A x+1 ~ is CB A x+1 - iscaA^B 2 + is BA A x - x C 2 

- isab{ Yl A p CA q AC) + is AC ( A p BA q AB) (28) 

p,q>o,p+q=x— 2 p,q>o,p+g=x— 2 
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where, since A, B and C satisfy the su(2) commutation rela- 
tions, we have that s BC = sign([£>,C]) = —s CBl Sba = 
sign([B,A}) = -s AB , and sac = sign([A, C] = -s C a- In 
each of the sums in (I28> . the various terms are of the same or- 
der (x + 1) and they differ from each other only for the order 
in which the A's and the B's(C's) appear. Since we want to 
apply the Casimir relation we have to first apply the ordering 
process discussed above to the terms in the sums to obtain: 

^2 A p BA q AB = nA x - l B 2 + p( x ) (29) 
I 



We now use the Casimir relation d26l to substitute (B 2 + C 2 ) 
I 



Since by hypotesis we already have all the terms of order < X 
we can use them to extract the desired term A x+1 . 

The next step of the procedure is based on the computation 
of commutators of the following form: 

[Ax-"i +1 ,Brt +1 ] = is A B A p B t CB s A\ 

P+9=X-7;*+s=7 

(33) 

where p,q,t,s > and 7 € {1, .., X — 1}. The terms in the 
sum have all the same sign sab and have order X + 1. Each 
term of the sum can be obtained by converting all the others 
into it by mean of an ordering process and then by subtracting 
the resulting linear combination of elements of order x- By 
changing 7 6 {l,..,x — 1} we can obtain all the terms in 
which the operator C appears just once; a prototype of these 
terms is 

A X ^B^C. (34) 

If we now start from these terms and repeatedly use the 
Casimir relation d26i we can get all the terms like 

A x+i-i-i B i c i (35) 

where I is an odd integer and now 7 £ {2, — 1}. In 
fact, suppose the desired term is A x ~ 2 ~ 7 B" 1 C 3 , then one can 
apply the Casimir relation to the sum of two appropriate terms 
of the form (I34> and obtain: 

A X ~ 2 ~^B^C(A 2 + B 2 ) = j(j + l)A x - 2 -" 1 B 1 C - 

_ A^^B^C 3 , (36) 



and 

53 A p CA q AC = nA^C 2 + p( X ) (30) 

p,q>0,p+q=x-2 

where n is the number of terms of each sum and p( X ) repre- 
sents a polynomial of order x- In this wa y e q- BH} can be 
rewritten as: 



(31) 



in l l3"H and we finally obtain: 



(32) 

I 

where j(j + l)yl x_2_ ' > B 1 C is of order X — 1 and can then 
be canceled, just as the terms of order x that would appear 
whenever an ordering process is applied to obtain the desired 
term. 

The next step of the procedure is based on the computa- 
tion of the commutator [A p , A a B l3 C] and on the use of the 
Casimir relation. By properly choosing the exponents p, a, /? 
the result of the commutator is the sum of terms of order x + 1 
and can be written in the following compact way: 

[A p ,A a B C] = imsABA^^B^^C 2 + 

+ ins A cA a+p - 1 B' 3+1 + p(x) (37) 

The term insAc^ a+p ^ 1 B l3+1 is the result of the ordering 
process on the n terms in which A appears a+p— 1 times and 
B appears (3+1 times. Analogously imsABA a+p ~ 1 B l3 ~ x C 2 
is the result of the ordering process on the m terms in which 
the operators A, B and C appear a+p—1, (3—1 and 2 times 
respectively. The term p(x) takes into account the terms of 
order X produced by these ordering processes. 

The procedure is based on the following steps. First we use 
the Casimir relation on A a+P ~ 1 BP+ 1 to get: 

A a+p~l B (3+l = A a +p -l B f3-l j( j + ^ _ 

- A^^B^A 2 - A a+p - x BP- l ClW) 

Then if we subsitute (I38i in J37i we can write: 

[A p ,A a B C] = i{m + n)sA B A a+p - 1 B' 3 - 1 C 2 + 

+ msAcA^^B^A 2 + p{ X ) (39) 



[A X ~ 1 B,CA] - {A X ~ 1 C,BA} = i2s B cA x+1 + i(n + l)s BA A x - 1 (B 2 + C 2 ) + p( X ) 

I 



[A^B, CA] - [A^C, BA] = i(3 + n)s BC A x+1 + i{n + ^sbaA^ 1 ^ + l)j + p( X ) 



8 



where we have taken into account that sac = —sab- We 
first note that for a + p — 1 = % — 1 an d P = 1> since we 
already have A x+1 , it is possible to generate all the terms of 
the type A X ~ Y C 2 and, by simply exchanging the role of C and 
B, A X ~ 1 B 2 . If we now set a + p — 1 = \ — 2 and (3 = 2, we 
can obtain the terms of the form A X B from d39i and A X ~ 2 B Z 
from i31i . In fact, in both cases all we have to do is to cancel 
out the terms like A X ~ 2 BC 2 that has been already obtained, 
seel34l 

For j3 > 3 the procedure use at each step the results ob- 
tained in the previous ones, that correspond to lower values of 
(3: first use (1391 to obtain the terms like 

Ax-PB^C 2 , (40) 

then use d37i to get the terms of the form 

A x ~ p B p+l . (41) 



The last part of the procedure is basically the same 
used to obtained (I35> : we start from the terms of the type 
AX-Pb^C 2 and, by repetedly using the Casimir relation 
we can get all the terms of the type 

A x+x ~ l -^WC l (42) 

where now I is an even integer. For example to obtain 
B~> 'C 4 we can write: 

Ax- 3 -~<B~<C 2 (A 2 + B 2 ) = A x - 3 - 7 B 7 C 2 j(i + 1) - 

- A x - 2 -^B^C 4 (43) 

where the first term on the right hand side is of order \ ~ 1 
and then can be canceled. 

If we now collect the results <32j , (f34j , 03 , <|40| , <f4ll> and 
(142 \ we see that we have succeded in giving an operative pro- 
cedure to build all the element of the enveloping algebra of 
su(2), starting from the set of available controllable interac- 
tions, i.e., 2^ = {X, Z, Z 2 }. 

VI. MASSLESS BOSON MODES 

Another interesting physical example is given by a set of L 
bosonic modes each representing the state space of massless 
bosonic particles. The quantum state-space in this case given 
by H.f = ®i = \hi, where h — {\n) } n= o,. .,oo- Independent 
modes of the electromagnetic field e.g, photons with different 
polarization, provide a prototypical and ubiquitous physical 
example of this kind of systems. 

Each of single-mode Hilbert spaces hi is an irreducible rep- 
resentation of the Weyl-Heisenberg algebra 

h(l) = {*,4,l}, [c h 4] = i,(i = l,...,L). (44) 

According to the general theory of section III, in order to 
achieve universal control over Tip one has to be able to gen- 
erate, as in theformer SU(2) case, the wole Hermitean part 



of the enveloping algebra of (g)f =1 hi(l) i.e., all the possible 
polynomials in the Cj's and cj's. Notice that this enveloping 
algebra is now given in an infinite-dinsional representation, 
whereas in the former SU (2) bosonic example the number su- 
perselection allowed us to focus on finite-dimensional irreps 
(the Fock sectors). 

It is a well known result in quantum control theory fl5ll that 
if is able to control -for example for the mode i- the following 
set of interactions 

l Ai ^{c^clic^^clf.n^inA 2 } (45) 

then one has universal control on the mode-system, i.e., one 
is able to generate, via commutation, linear combinations of 
Hermitian polynomials in Ci and c\ of any degree. The pres- 
ence of the non-linear tern (n^) 2 is essential, because it en- 
tails to increase via commuation the order of the polynomials 
already available. 

Furthermore, if in addition to I Ai one has the control over 
the set I a and of two-modes interactions: 

' \ i = {'v. -'V-!' 

then one has the universal control over the whole two-modes 
system,i.e. one can generate, via commuation, linear combi- 
nations of Hermitian polynomials in Cj,c|, Cj and cj of any 
degree 1 15]. 

What we first want to show here is that it is sufficient to 
control the non-linear terms of only one of the bosonic modes, 
say (ni) 2 , in order to have universal control over each of the 
modes i £ {2, .., L}. If we suppose to be able to control the 
following set of interactions: 

l Al ={c 1 ,c t 1 ,(c 1 ) 2 ,(c t 1 ) 2 ,n 1 ,(n 1 ) 2 }, (46) 



x Ai = {a,cl (c,) 2 , (ct) 2 , m}, i e {2, .., l}, (47) 

and 

then we see that in order to reach our goal we have to prove 
that we are able to generate each of the terms nf, i £ 
{2, .., L}. This can be done by using the following commua- 
tors: 

[[nl,T u ],T u ] = ~8n ini + 2c 2 (cJ) 2 + 2c 2 (c t 1 ) 2 + P; (49) 
[niJni.Knf.Tiil.Tn]]] = 8c 2 (4) 2 + 8c 2 (4) 2 ; (50) 

[[mriijTu^Tu] =n% + Q, (51) 

where Tu = c\ci + c\c\ and P, Q are polynomials in terms 
already available. First the procedure requires the use of d50i 
in order to extract from ( I49> the term n\fii. Then the latter 
is used in ( 15 1> to generate the desired term nf. We have thus 
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shown that it is possible to "propagate" the non linearity n\ to 
all the other modes. 

The control over the sets of interactions J45L ( 147 \ . J48t 
and of the terms rif allow 1 15] universal control over any of 
the two-modes subspaces hi (8 hj. In fact with the previous 
hypotesis, given a two-modes operator Oi ® one can al- 
ways generate all the possible operators Oi ® Pj, where Pj 
is a generic linear combination of Hermitian polynomials in 
Cj and Cj of any degree; this can be done by judiciusly com- 
muting Oi (g> Aj with operators 1 ® acting on the j'th 
mode only. Analogous arguments can be used to generate 
multi-modes operators, and this means that the sets of inter- 
actions ( I45> . ( I47K ( 1481 allow universal control over the whole 
L-modes state space Tip- 

VII. AUXILIARY BOSON 

In this section we will discuss a technique for implement- 
ing effective Hamiltonians in a group-algebra C/C which in- 
volves the use of an auxiliary bosonic mode. This an ap- 
plication of thegeometric -phase operators discussed in Refs 

asms 

We consider a quantum state-space Ti which is an irrep 
space of the group JC. The representatives of the group ele- 
ments have the form g = cxp(iH g ) E End(7i) (H g = H^). 
Our aim is to generate an effective dynamics governed by an 
Hamiltonian in uC/C, of the form g + g^ or i(g — g''). These 
two operators, in terms of the Hamiltonians H g , are clearly 
proportional to cos(H g ) and sin(_ff g ) respectively. Now the 
point is that it is known how to enact evolutions associated to 
this kind of highly non-linear operators by resorting to a tun- 
able coupling with an ancillary bosonic mode. This coupling 
has to be of the form K g :— 9H g ® n where n — a^a is the 
boson number operator ([a, a^] = 1) Moreover one has to be 
able to control the bosonic mode itself in order to enact the 
dynamics described by the displacement operators 

D(a) := exp(aa - aa'), (a € C). (52) 

When one has at disposal such ingredients the following rela- 
tion shows how to generate the desired effective Hamiltonian 
over Ti while leaving the bosonic mode unaffected 

D{-(3)D(-ae l9H <>)D(!3)D(ae ieH °) = e ^M^(<>H s +4>) 

(53) 

Where_ D{ae l8H = ) = e l6nH <> D(a) e~ l6H <> and := 
arg(a/3). 

For a finite order group JC the most general element in uC/C 
has the form J2 g (^g C0S (H g )+n g sin(iJ g )), (A g ,/i g S H).It 
follows that if one has available all the couplings K g the simu- 
lating sequence (I53> allows to generate all possibile Hamilto- 
nians in the group algebra of JC. Even when just some, possi- 
bly just one, K g coupling is actually available on may still be 
able to generate the whole uC/C, by generating commutators 
with other simple Hamiltonians at disposal. 

In order to exemplify this latter point let us consider the 
continous group SU(2) and its irreps corresponding to total 



angular momentum J. Since the irrep has dimension 2 J + 1 
one has e.g., that sin(6S z ) = jfkLo Fj(0){S z ) k where the 
Fj are real- valued functions. The implementation of the 
trascendental function sin in finite-dimensional spaces just 
amounts to generating a high order polynomial. From the pre- 
vious analysis we know that this, along with the capability to 
switch on and off simpler Hamiltonians e.g., S X ,S V , may be 
sufficient for generating all the Hamiltonians in C/C. 



VIII. CONCLUSIONS 

In this paper we thoroughly analyzed the dynamical- 
algebraic approach to universal quantum control introduced 
in II Oil. The basic idea is to decompose the quantum state- 
space TL encoding information into irreducible sectors and 
subsystems associated to the group of available evolutions. 
When this group coincides with the unitary part of the group- 
algebra C/C of some (non-abelian) group JC then universal 
control is achievable over the /C-irreducible components of 
H. Physically this can be done by turning on and off a 
generic pair of Hamiltonians in C/C. We then applied this 
general strategy to different kind of bosonic systems. We 
first considered massive bosons in a double-well and showed 
how to get universal control over all the finite-dimensional 
the Fock sectors by the control of single-well frequency, tun- 
neling and a non-linear term. We then discussed the multi- 
mode massless case showing again that linear (in the optics 
sense) Hamiltoniain along with a single non linearitity allows 
one to generate the whole (infinite-dimensional) multi-mode 
Heisenberg-Weyl enveloping-algebra and so to achieve uni- 
versality. Though some similar results have been already de- 
rived in other situations it is remarkable the fact that they pro- 
vide a conceptually straightforward and enlighting exemplifi- 
cation of the general control-theoretic framework introduced 
in Ref. IT Oil . Notice that in those bosonic examples we gave 
a constructive way to generate effective interactions, out of a 
simpler ones, by means of higher order commutators. This 
in turns provides a explicit simulation sequences for complex 
interactions e.g., high-order non-linearities in terms of read- 
ily available ones. Finally we discussed how to use an aux- 
iliary bosonic mode coupled with finite-dimensional systems 
to generate the high-order non-linearities necessary to realize 
universal control. We believe that the results presented here 
shed further light over the strengh and reach of the universal 
quantum control approach introduced in 1 10]. 
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